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ABSTRACT 



We present a detailed stellar population analysis of 27 massive elliptical galaxies 
within 4 very rich clusters at redshift z ~ 0.2: A115, A655, A963 and A2111. Using the 
new, high-resolution stellar populations models developed in our group, we obtained 
accurate estimates of the mean luminosity-weighted ages and relative abundances of 
CN, Mg and Fe. We derived the age, [CN/H], [Mg/H], [Fe/H] and the abundance ratios 
[CN/Fe] and [Mg/Fe] as functions of the galaxy velocity dispersion, a. We have found 
that [CN/H] and [Mg/H] are correlated with a while [Fe/H] and Log(age) are not. In 
addition, both abundance ratios [CN/Fe] and [Mg/Fe] increase with a. Furthermore, 
the [CN/H]— a and [CN/Fe]— a slopes are steeper for galaxies in very rich clusters than 
those in the less dense Virgo and Coma clusters. On the other hand, [Mg/H]— a and 
[Mg/Fe]— a slopes keep constant as functions of the environment. Our results are com- 
patible with a scenario in which the stellar populations of massive elliptical galaxies, 
independently of their environment and mass, had formation timescales shorter than 
~ 1 Gyr. This result implies that massive elliptical galaxies have evolved passively 
since, at least, as long ago as z ~ 2. For a given galaxy mass the duration of star for- 
mation is shorter in those galaxies belonging to more dense environments; whereas the 
mass-metallicity relation appears to be also a function of the cluster properties: the 
denser the environment is, the steeper are the correlations. Finally, we show that the 
abundance ratios [CN/Fe] and [Mg/Fe] are the key "chemical clocks" to infer the star 
formation history timescales in ellipticals. In particular, [Mg/Fe] provides an upper 
limit for those formation timescales, while [CN/Fe] apperars to be the most suitable 
parameter to resolve them in elliptical galaxies with a < 300 km s _1 . 

Key words: Galaxies: abundances - Galaxies: clusters: general - Galaxies: evolution 
- Galaxies: formation - Galaxies: stellar content 



1 INTRODUCTION 

Stellar Population analysis has recently shown increasing 
power in constraining galaxy formation scenarios. This 
follows the relative success of CDM cosmologies, giving 
rise to the hierarchical build-up of g alaxies as a con- 
sequence of the initia l cond i tions (e.g . Press fe Schechtei] 



19741: IWhite fe Frenkl Il99ll ; ISomerville fe Primackl 119991 : 
De Lucia et alj|2006h . Many structural and dynamical prop- 
erties of galaxies in clusters are explained in this scenario, 
though problems remain: for instance, the absence of the 
pred icted mass cusps in the centres of ellipticals and bulges 
fe.g. IValenzuela et al.ll2005f ). and the prediction of far more 
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satellite galaxies than those observed l|Klypin et al.l 1 19991 : 
iMoore et aTlll999l ). 

Stellar populations offer a fossil record of the star for- 
mation and chemical evolution of galaxies, most clearly in 
elliptical galaxies, thought to have formed by the merging of 
discs, so stellar population studies provide very strong con- 
straints on these galaxy formation scenarios. For example, 
it is particularly hard to reconcile the hierarchical models 
with the result that massive galaxies show significantly older 
mean luminosity weighted ages than their smaller counter- 
parts (|Kauffmann et alj|200a lYamada et alj|2006h . 

A basis for population analy sis is the colour-magnitude 
relation (CMR). ISandagel (1 19721 ) showed that luminous el- 
liptical galaxies are redder than fainter ones. The mono- 
lithic explanation is that more massive galaxies produce 
more metals at an early phase (jKodama fe Arimotol [l997l : 
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IChiosi fc Carraroll200l lMatteuccill2003T l. implying old stel- 
lar populations in elliptical galaxies. On the contrary, in 
the hierarchical explanation the metals are produced over 
an extended period in merger induced star formation, 
so an intermed i ate a ge population should be detectable 
IIDe Lucia et al ] l2006h . The results so far a e ambiguous: 



Terlevich et al, Q 19991 ) found evidence for the former, while 



J0rgensenl l| 19991 ') and lTrager et all l|2000l) for the latter 



It is hard to separate age from metallicity effects us- 
ing the CMR (|Worthev et all 1 19941 ) and this bedevils all 
analyses of unresolved stellar populations. Many studies of 
clusters have used th e pioneering Lick/IDS spectral indices 
l|Worthev et al.lll994l ) whose reliablity is limited by their res- 
olution dependence as uncertain corr ections for broadening 
and instrumental effects are needed IjWorthev fc Ottavianil 
1997). One can make some advance by using modelled in- 
tegrated spectra, which can be broadened to match the 
galaxy velocity dispersion, sigma (although indices are still 
used in this method, one has more flexibility that with the 
fitting functions required in the earlier technique). Mod- 
els employing the use of integrated spec tra at re l atively 
high resolution have been developed by IVazdekisl (|l999l ) 
which allow one to effect a better se paration of age and 
metallicity than with the Lick indices ( Vazdekis fc Arimotq 
Il999l ; IVazdekis et~ai]|200ll ; iKuntschner et al.ll2002l ). Another 
advantage we have in the present work is that we have 
been able to employ the new stellar population synthe- 
sis models of Vazdekis et al. (2006, in preparation) which 
are improved versions of the earlier IVazdekisl l|l999l ) mod- 
els, based on a greatly extended stell ar spectral library 
, MILES (ISanchez-Blazquez et all 12006? ) . The new library 
contains spectra of over 1000 stars, convering unequalled 
ranges in stellar metallicity and spectral class in the wave- 
length range lambdalambda AA 3500 — 7500 A, at a resolu- 
tion of 2.3 A (FWHM). 

This new generation of models has shown considerable 
ability to break the age metallicity degeneracy and has been 
succesfully applied to a number of elliptical galaxies se- 
lected using the CMR of the Virgo cluster. Once the age 
metallicity degeneracy is resolved, it is fairly easy to derive 
the abundanc es of individual eleme nts using these models. 
For example, [Vazdekis et al.l l|200ll ) have found that lumi- 
nous and moderately-luminous ellipticals (Mv< — 18.3 mag; 
a ^ 100 km s _1 ) are essentially old (t > 8-10 Gyr) without 
a clear trend with luminosity. [Mg/H] correlates tightly with 
velocity dispersion, while [Fe/H] is showing only a modest 
trend, implying that the CMR is driven essentially by metal- 
licity (in particular, by the abundance of the a-enhanced 
elements) and that star formation in massive ellipticals vir- 
tually had stopped before SNIa started to contaminate the 
ISM with iron-group elements. In short, star formation was 
intense and short-lived (^ 1 Gyr). Admittedly, however, the 
number of ellipticals observed in this way is small and the 
results may not be sta t istical ly significant. Contrary to these 
findings, llrager et all (|2005T ) found a few young m oderately 
massive ellipticals in Coma. ICaldwell et al.l (|2003h analyzed 
integrated spectra of 175 nearby early-type galaxies (mix- 
ture of field and Virgo E, dE, SO, and dSO's) to find a trend 
between the population age and the velocity dispersion, such 
that low-er galaxies have younger luminosity-weighted mean 
ages. It is therefore not yet clear if cluster ellipticals with 



a ^ 100 km s _1 are uniformly ol d or are progressivel y 
younger towards less massive objects l|Bernardi et aljfeood ). 

One of the key pieces of information is to see how abun- 
dance ratios in early-type galaxies vary with the cluster 
masses (as sampled via their richness class or their X-ray 
luminosities, for example) and, in this way, obtain vital in- 
sights into galaxy assembly timescales, as functions of t heir 
environment (e.g. IWorthevlll998l ; ICarretero et alj|2004 ). In 
particular, overabundances of [Mg/Fe] compared with the 
solar ratio have been found in ma s sive elliptical galaxie s 
jPeletierl Il989l ; IWorthev et ail 1 19921 ; IVazdekis et ailll997f l. 
These have been interpreted via several possible scenarios 
based on th e fact that Mg is mainly produce d in Type II 
supernovae ([Faber et al ] |1992| ; lMatteucci|[l994l ) , and include 
different star formation rates (SFR) and a time dependent 
initial mass function (IMF). 

Bas ed on SPSS data, we followed this approach pre- 
viously (ICarretero et alj [2004?). testing the evolution of the 
ratios between CN, Mg and Fe as a function of the cluster 
X-ray luminosity. We tentatively concluded that early-type 
galaxies are fully assembled on timescales around the epoch 
of the massive release of CN into the ISM. This timescale is 
shorter than predicted by c lassical hierarchical models (e.g. 
iKauffmann fc Charlotlll998l ). and thus it is vital to test this 
conclusion. 

In this paper, we concentrate on the stellar popula- 
tion analysis of 27 elliptical galaxies distributed in 4 rich, 
intermediate-redshift (z ~ 0.2) Abell clusters: A115, A655, 
A963 and A2111. Selected cluster galaxies have been mor- 
phologically and photometrically studied in the literature, so 
we focused on those galaxies already classified as ellipticals. 
We determined mean luminosity-weighted ages and a num- 
ber of abundance ratios of species such as Fe, CN and Mg; 
and their relation with the velocity dispersion of the galax- 
ies. Finally, we compared our results with those obtained by 
other authors in galaxies belonging to less rich clusters such 
as Virgo and Coma. 

The outline of the paper is as follows: Section 2 presents 
the observations, sample selection and data reduction. Sec- 
tion 3 shows the relations found between the abundances of 
selected elements and the velocity dispersion, and the com- 
parison with results obtained in less dense environments. 
Section 4 presents the discussion of our results. Finally, Sec- 
tion 5 summarizes our conclusions. 

Where necessary, we adopt a cosmological model with 
H =71 km s" 1 Mpc" 1 , f2 m =0.3 and f2 A =0.7. 



2 OBSERVATIONS AND DATA 
2.1 Sample selection 

We have studied a set of 27 elliptical galaxies belonging 
to 4 very rich (ric hness class 3 and Lx ~ 10 45 erg s -1 ; 
lEbeling et all 1 1998! ) Abell clusters: A115, A655, A963 and 
A2111. Table [1] presents detailed information on each clus- 
ter. We have selected these clusters because they represent 
very dense environments, since our aim is to compare the re- 
sults of the stellar population analysis performed on galax- 
ies in very rich clusters with the results already obtained 
in less rich clusters, su ch as Virgo (richn ess class 1 and 
L x = 0.3 x 10 44 erg s- MLedlow et al.ll2003l ) or Coma (rich- 
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Figure 1. Field of Abell 655. We have superimposed the MOS 
mask #1 used during our observations of this cluster. 



ness class 2 and L x = 1.8 x 10 44 erg s'^ lLedlow et alj|2003h 
clusters. 

A115 presents 2 main X-ray substructures, which sug- 
gests that this cluster is not relaxed but in an inter- 
mediate dynamical stage of merging prior to virializa tion 
jBeers et al.lll983l ; [Gutierrez fc Lopez-Corredoiral l2005j l A 
similar scenario can be traced for A655, w hich presents var- 



ious substructures ( Strazzullo etal 



20051). and for A2 111. 



which shows a double-peaked core jjeltema et al.ll2005h . A 
very different situation is found for A963 which, even though 
a high Lx cluster, is a quite relaxed system with some pos- 
sible very small-scale structure in the core (|Jeltema et al.l 
120051 ) . 

The selected galaxies are the brightest galaxies in their 
respective clusters (15.1 < rav < 17.2) and have been al- 
ready mor phologically analyzed in the literatur e (see Table[T| 



for details; iFlin et alj|l995t iFasano et ail l200d ; iRakos et al. 
l200d ; iBalogh et alJ^OoC Thanks to this previous work we 



were able to focus our observations on those galaxies already 
classified as ellipticals. The velocity dispersion values of the 
galaxies in our sample cover the range 220 km s _1 < a < 
450 km s _1 . 



2.2 Observations 

Multi-slit spectroscopy was carried out with the DOLORES 
spectrograph on the 3.5m Telescopio Nazionale Galileo at 
the Roque de los Muchachos Observatory (La Palma, Spain) 
in 3 runs. Table [1] details each observing run. DOLORES 
CCD scale is 0.275 arcsec pix -1 which yields a field of view 
of about 9.4 x 9.4 arcmin. We designed Multi Object Spec- 
troscopy (MOS) masks containing ~ 6 slits each. As an ex- 
ample, Figure [1] shows the field of A655. We have overim- 
posed the image of one of the 2 MOS masks used to obtain 
the spectra from the galaxies of this cluster. 

For all masks, in all observing runs, we used the same 
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Figure 2. Top: An example of one of our spectra. The black line 
corresponds to the observed spectrum of a representative galaxy 
of A115 (<r = 330 km s" 1 ). The red line corresponds to a model 
of Vazdekis et al. (2006, in preparation) with t = 15 Gyr and 
solar metallicity. Bottom: Residual after subtracting the model 
from the observed spectrum. 



instrumental setup: the medium-resolution MR-B#2 grism, 
in the wavelength range AA S3500-7000 A, and 1.1" width 
slits. The dispersion was 1.7 A pix with an instrumental 
resolution of 8 A (FWHM). 

We performed exposures of 1800 seconds, interspersed 
with He+Ar arc-lamp spectra in order to correct for possi- 
ble changes in the wavelength calibration during the night. 
Table [1] presents the total exposure time for each mask. In 
addition, we observed a number of G- and K-type stars to 
be used as templates for velo city d ispers ion measurements. 
Finally, flux standards from were observed to 

correct the continuum shape of the spectra. 



2.3 Data Reduction 

The standard data reduction procedures (flat-fielding, wave- 
length calibration, sky subtraction and fluxing) were per- 
formed with I RAF packages. In addition, we used REDmE 
l|Cardiel| [l999) to remove cosmic rays, since this reduction 
package allows full control of this process. 

Once reduced, we coadded our 2D spectra in the spa- 
tial direction to obtain an average value of S/N~ 40 in 
the full sample. Considering the physical distance to each 
cluster, we integrated the light of the galaxies out to a 
galactocentric radius comparable to th e galaxies already an- 
alyzed in Virgo and Coma c lusters (|Vazdekis et al.l |2004| ; 
ISanchez-Blazquez et al.l l2006). In particular, we considered 
an equivalent aperture of 5" at a redshift z = 0.016. This 
way, we were able to establish direct comparisons between 
the our spectral index measurements and those performed 
in these clusters. 

Figure [2] shows an example of the final spectrum of a 
galaxy in A115 with a — 330 km s _1 . We have overplot- 
ted a synthetic spectrum of a model of age 15 Gyr and 
solar metallicity. The bottom panel presents the residuals 
obtained when subtracting the model spectrum from the 
observed one. 



3 ANALYSIS AND RESULTS 

Once the reduced spectra for all the galaxies in our sample 
were obtained, we proceeded to measure their velocity dis- 
persion and different spectral indices, which are the basis of 
our stellar population analysis. 
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Tabl e 1 ■ Logbook of the o bservations and main properties of selected Abell clusters. X-ray luminosity was obtained 
from lEbeling et al.l <j 199Sf ) ■ N ga i column contains the number of galaxies studied for each c luster. " Ref . " colu mn 
indicates the refe ren ce from which we ha ve obtained the morpho logy of analysed galaxies: l: lRakos et al. | (200Cft ; 2: 
iFlin et alj <1995ft : 3: lBalogh etall d2002Tl ; 4: lFasano etHI i200Cl) . 
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a iexp per mask (2 in total). 

k The observing program was interrupted for several hours during two nights due to the occurrence of 2 GRB 
warnings: GRB 050505 and GRB 050508. 



3.1 Velocity Dispersion Measurements 

To measure the velocity dispe rsion for each galax y we fol- 
lowed the method described in lDavies et al] l| 1993t h using a 
series of G- and K-type star spectra obtained with the same 
instrumental configuration as the galaxy spectra. 

The spectrum of a galaxy is broadened due to: 1) the 
intrinsic velocity dispersion of the galaxy stars, a, and 2) 
the instrumental resolution, ai nstr . Thus, the final observed 
velocity dispersion, <7 b s , is given by: 



2 _ 2 , 2 
G obs — ° + f instr 



(1) 



On the other hand, a stellar spectrum is affected only 
by the instrumental resolution, since a single object has 
zero intrinsic velocity dispersion. In this case, a = and 
o"obs = Omstr- This method of measuring the intrinsic ve- 
locity dispersion of a galaxy entails broadening the stellar 
spectrum using a series of Gaussian niters with increasing 
standard deviation, a. 

Following this procedure, we estimated the intrinsic ve- 
locity dispersion for the full set of galaxies. Table [2]presents 
the measured a values. Note that our sample includes mas- 
sive and very massive objects since we are covering the range 
220 km s _1 < a < 450 km s -1 . 



3.2 Ages and Metallicities 

To derive mean luminosity-weighted ages and metallicities, 
we compared selected absorption line strengths with those 
predicted by the model of Vazdekis et al. (2006, in prepa- 
ration). This model provides flux-calibrated spectra in the 
wavelength range AA 3500 - 7500 A, at a resolution of 2 A 
(FWHM) for si ngle-burst stella r populations. This model is 
an extension of IVazdekisI J 19991). employing the new empir i- 
cal stellar spectral library of lSanchez-Blazquez et al] (2006). 
This way, we can transform synthetic spectra to the reso- 
lution and dispersion of the galaxy spectra instead of the 
opposite, as required when working in the Lick/IDS system. 

Once the model spectra are transformed to yield the in- 
trumental conditions of resolution and dispersion of the ob- 
served spectrum, we measure pairs of indices in both sets of 
data (observed and synthetic). These pairs are formed by an 
index which is sensitive almost only to variations in age, and 
other which is sensitive almost only to variations in the abun- 



dance of a given element. Then, we obtain nearly orthogo- 
nal grids of age versus metallicity which can be compared 
directly with the values of the same pair of indices measured 
in the observed spectrum. Figure [3] illustrates this method. 
The vertical axis shows the age indicator H/3 while the three 
horizontal axes show the metallicity indicators Fe4383, CN2 
and Mg2, respectively. By comparing the model grids with 
the measurements in the observed galaxy spectrum we can 
estimate metallicities and mean luminosity-weighted ages. 
Note that the full set of synthetic model spectra is broad- 
ened to the velocity dispersion value of each galaxy, i.e. we 
create individual grids suitable for each galaxy. This way we 
do not need to correct for the velocity dispersion effects, to 
match a fixed resolution as is the case of the commonly used 
Lick/IDS analysis. 

The selec ted metallicity indice s were CN2, Mg2 and 
Fe4383 (c. f. IWorthev et all 1 19941) . We used H/3 as the 
age index in all cases. We chose these features because 
of their low sensi tivity to variations in signal-to-noise 
l|Cardiel et al] [20031 ) and velocity dispersion (we estimated 
A(index)/ index < 0.15, for A(cr) = 300 km s _1 ). This way 
we minimize possible radial variations in the index value, as 
a may vary as a fun ction of the galactocentric radius (see 
ICarretero et al.ll2004l ). 

We will refer to the metallicities derived in the diagrams 
CN2-H/3, Mg 2 -H/3 and Fe4383-H/3 as [CN/H], [Mg/H] and 
[Fe/H], respectively. Since the CN2 index is strongly domi- 
nated by C and N, the Mg2 ind ex is governed by Mg and the 
Fe4383 index is driven by Fe (|Tripicco fc Bellill995l ). these 
metallicities must be close to the real relative abundance of 
each element. Given the values of the former relative abun- 
dances, we can obtain estimates of the abundance ratios 
[CN/Fe] and [Mg/Fe] for each galaxy. 

The bulk of our galaxy sample lies with in the bottom- 
right corner of the model grids shown on Figure [3] However, 
note that an extrapolation of the model grids is required 
for some galaxies to obtain the abundances of CN and Mg, 
since the models extend only to [M/H] = 0.2. Also, it is 
worth noting that the data points for certain galaxies fall 
below the model grids, which can be attributed to the fact 
that absolu te age determination is subject to model zero- 
points (see I Vazdekis et al.l I2OO1I ; ISchiavon et all 120021 ) . We 
neglect the possible effect of nebular emission on H^ since 
no [O III] A5007 emission is detected in our galaxy spectra. 
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Figure 3. Examples of the method used to derive ages and metallicities. The plot shows the age index H/3 versus the metallicity 
indicators Fe4383, CN2 and Mg2 for the grid of models of Vazdekis et al. (2006, in preparation). Horizontal dotted lines represent models 
of constant age whose ages are labelled in the plot. Vertical dashed, solid, dotted, dashed-dotted, dashed-dottcd-dotted and long-dashed 
lines represent contours of constant [M/H] = +0.2, 0.0, —0.4, —0.7, —1.2 and —1.7 dex, respectively. Diamonds, open circles and squares 
indicate the values for a representative galaxy of A115, A655 and A963, respectively, with a common a = 310 km s _1 . Note that all 
synthetic model spectra were broadened to the instrumental resolution of the observed spectra and the velocity dispersion of these 
particular galaxies (<r obs = 350 km s _1 ). 



Neverthe less, assuming an upper emission correction of 
~ 0.5 A (|Kuntschner et al]|2002l ). although this would give 
rise to a significant reduction in the mean age of the stel- 
lar populations of the oldest galaxies, the net effect on the 
abundance ratios would be no more than ~ 0.05 dex, for the 
most affected galaxies in a cluster. This correction for the 
H/3 index is larger than that obtained by varying the model 
prescriptions, e.g. with q-enhanced iso chrones and atomic 
diffusion included (Vazd ekis et alJfeOQll ). used to derive ages 
in better agreement with the current age of the Universe, for 
the oldest galaxies in our sample. 

Finally, Table [5] lists the estimates of the values of the 
relative abundances [CN/H], [Mg/H] and [Fe/H] and the 
mean luminosity-weighted ages for each galaxy in our sam- 
ple. 

3.3 Results 

3.3.1 Correlations with the Galaxy Mass 

Once we have obtained estimates of the relative abundances 
of selected elements, and the mean luminosity-weighted age 
for each galaxy, we studied their relation to the galactic 
mass, via the velocity dispersion. Figure [3] presents the re- 
lations [CN/H]-cr, [Mg/H]-cr, [Fe/H]-cr and Log(age)-<r 
for our full sample. The solid line in each panel corresponds 
to an error-weighted linear fit to the data. We can see that 
the abundances [CN/H] and [Mg/H] increase with a. On the 
contrary, we find no clear trends with a for [Fe/H] or the 
age. 

In order to compare our results from galaxies in very 
rich clusters with those from galaxies within less dense en- 
vironments, we have also plotted on Fig ure U the results 
obtained bv lSanchez-Blazquez et all (|200€ ). who analyzed a 
sample of 98 early-type galaxies belonging to the Virgo and 
Coma clusters. The dotted line in each panel corresponds to 
the relation found for the Virgo cluster galaxies, while the 
dashed line corresponds to the Coma cluster galaxies. Note 
that the sigma coverage is limited to a < 350 km s _1 for 
the Virgo and Coma cluster samples. 

In agreement with our findings, [CN/H] and [Mg/H] are 



correlated with a for the Virgo and Coma cluster galaxies. 
Also, no clear trend is found for [Fe/H] —a. When we con- 
sider the age, the result differs from Virgo to Coma: while 
Virgo galaxies show a clear trend with a, none is found for 
the denser Coma cluster. In particular, we observe that the 
behaviour of the Log(age) — a relation is very similar for the 
Coma galaxies and for our sample. 

Note that, in the bottom panel of Figure 0] there is an 
offset between the Log(age)— a relation for our cluster sam- 
ple and the Virgo and Coma ones. This is not a critical issue 
since we are only interested in the slopes of the correlations, 
and not in their absolute values. Nevertheless, the analysis 
of the relative age values does give useful information. This 
is also the case for derived metallicity values: we do not ex- 
pect to obtain absolute values of [M/H], since we based our 
analysis on the relative ones. 

We have focused so far on the relative abundances of 
selected elements. However, in order to obtain information 
concerning characteristic formation timescales for the bulk 
of their stellar populations, we computed the abundance ra- 
tios [CN/Fe] and [Mg/Fe]. Since each element was ejected 
into the ISM by different types of stars at different epochs, 
their relative abundance values will be very useful for con- 
straining the star formation history. Figure [5] presents the 
relations of the abundance ratios [CN/Fe] and [Mg/Fe] with 
the velocity dispersion in our sample. The solid line in each 
panel corresponds to an error- weighted linear fit to the data. 
We observe positive trends with a for both abundance ratios 
[CN/Fe] and [Mg/Fe]. 

As in the previous relations, we compared our results 
with those obtained in Virgo and Coma clusters. Figure 
[5] shows the results. The dotted line in each panel cor- 
responds to the relation found for Virgo cluster galaxies, 
while the dashed line corresponds to Coma cluster galax- 
ies. In this case, we also see correlations of both [CN/Fe] 
and [Mg/Fe] with a. We have also compared our result with 
the [Mg/Fe]— a relati on for Virgo cluster galaxies found by 
IVazdekis et al.l l|2004 ). which is represented in bottom panel 
of Figure [5] by a dashed-dotted line. 
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Table 2. Determined values of the velocity dispersion, mean luminosity-weighted age and relative abundances for our 
galaxy sample. 
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0.26 


0.76 


0.08 


0.66 


0.13 


0.18 


0.11 


A115 — 4 
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10.16 


0.29 


0.34 


0.11 


0.34 


0.13 


-0.21 


0.14 
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370 


10.01 


0.31 


1.23 


0.26 


0.63 


0.28 


0.08 


0.18 


A115 — 6 


330 


10.18 


0.22 


0.04 


0.12 


-0.09 


0.12 


0.03 


0.12 


A655 — 1 
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10.05 


0.36 


0.19 


0.29 


0.70 


0.22 


0.06 


0.33 
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0.17 
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0.10 


0.25 


0.17 
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0.47 
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0.31 
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-0.17 
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0.11 
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10.10 
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-0.69 
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A963 - 1 
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A963 - 2 
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0.87 
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A2111 - 3 
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0.64 
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10.19 
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1.30 


0.08 


1.30 


0.08 
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0.12 


A2111 - 5 


230 


9.85 


0.38 
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0.54 


0.35 


0.13 


0.11 


0.15 



3.3.2 Correlations with the Environment 

Figure [6] presents the slopes of the above Z— a relations 
(where Z represents any of the abundance ratios above) 
as a function of the cluster X-ray lum inosity, which is a 
Rood indicator of the cluster mass (e.g. lEttori et al.ll200j ; 
iDonahue et~al1l2003l ; Ishimizu et al.ll2003h and a useful rep- 
resentative of its richness class. This way we can compare 
the strength of the full set of relations as a function of the 
environment. Figures [6^, to[6}l present the slopes of the re- 
lations [CN/H]-cr, [Mg/H]-cr, [Fe/H]-(T and Log(age)-cr 
as functions of Lx, for our sample and for the Virgo and 
Coma cluster galaxies. For the slope of the [CN/H]— a rela- 
tion, we observe a clear trend of the slopes to increase with 
the density of the environment. On the contrary, when con- 
sidering the slopes of the [Mg/H]— a relation we do not find 
a clear trend with the environment, all slopes being almost 
constant and positive. In addition, the [Fe/H]— a relation 
presents slopes close to zero independently of the proper- 
ties of the host cluster. Different behaviour is found for the 
relation of the slope of Log(age)— a as a function of the en- 
vironment. In this case, our results suggest that this slope 
is anticorrelated with the X-ray luminosity. 

Finally, Figures [6^ and present the dependence of 
the slopes of the [CN/Fe]— a and [Mg/Fe]— a relations, on 
the cluster Lx- We can observe that, for both abundance 
ratios, the slopes are positive. In addition, the slope of the 
[CN/Fe]— a relation increases with the cluster mass, while 
the slope of the [Mg/Fe]— a relation is nearly constant, 
though positive, within the errors. 



4 DISCUSSION 

In this Section we will try to explain within a common sce- 
nario all the trends found. In particular, the correlations can 
be interpreted in terms of the different formation timescales 
for each chemical element, and the different star formation 
histories of elliptical galaxies as a function of their environ- 
ment. 

Magnesium is ejected into the interstellar medium 
(ISM) by Ty pe II supernovae ( S Ne II) on short timescales 
(< 10 Mvr; iFaber et all 11993 ; iMatteuccil 1 1994 ) . On the 
other hand, the iron-peak elements are the products of 
SNe la, which occur on timescales of ~ 1 Gyr. Be- 
tween the two extremes, although there are recent sug- 
gestions that most of the C comes from massive stars 
l|Akerman et alj|2004l ). C and N are mainly ejected into the 
ISM by low- and intermedi ate-mass stars (R enzim fc Volil 
Il98ll ; IChiappini et al.l l2003h . leading to CN formation on 
timescales longer than tha t for Mg but shorter than that 
for F e. Several authors (e.g. lEllis et aill 19971 ; IStanford et al.l 
1998) argue that early-type galaxies are old and passively 
evolving systems. The luminosity-weighted ages derived 
from our model grids confirm that the galaxies are signif- 
icantly older than the formation timescales of the different 
species. So, if we find substantial differences in the abun- 
dance ratios of these elements which depend on the physical 
properties of the environment, these must be due to the fact 
that galaxies are assembled on different timescales as a func- 
tion of their environment. 

We will focus our discussion on the following three as- 
pects: 
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Figure 4. Relations between the relative abundances [CN/H], [Mg/H], [Fe/H], and mean luminosity- weighted ages and the velocity 
dispersion of the galaxies, a. Diamonds: A115 galaxies. Open circles: A655 galaxies. Open squares: A963 galaxies. Open triangles: A2111 
galaxies. The solid line corresponds to the error- weighted linear fitting to our data. Dotted line and dashed line corresponds to the 
relations obtained for Virgo cluster galaxies and Coma cluster galaxies, respectively (taken from Sanchcz-Blazquez ct al. 2006). It is 
noteworthy that the slope of the relations increases with the environment density for [CN/H]— a, while it is almost constant in the case 
of [Mg/H]-<7 and [Fe/H]-o\ 



1) The comparison of the relative abundances and abun- 
dance ratios measured in galaxies of the same mass (i.e. 
equal velocity disperion) in different environments. 

2) The comparison of the stellar population properties 
measured in galaxies of different mass within the same en- 
vironment. 



3) The dependence of above relations as functions of the 
galaxy mass, with the environment properties. 

First, when inspecting Figures [4] and [5] we observe that, 
for a given a, the relative abundance [CN/H] increases with 
the cluster X-ray luminosity, whereas [Mg/H] and [Fe/H] do 
not. Something similar happens with the abundance ratios: 
while [CN/Fe] increases as a function of the environment 
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Figure 5. Abundance ratios [CN/Fe] and [Mg/Fe] versus the velocity dispersion of the galaxies. Diamonds: A115 galaxies. Open circles: 
A655 galaxies. Open squares: A963 galaxies. Open triangles: A2111 galaxies. The solid line corresponds to the error-weighted linear 
fitting to our data. Dotted line a nd dashed line corresponds t o the relations obtained for Virgo cluster galaxies and Coma cluster 
galaxies, respectively (taken from Sanchcz-Blazaucz ct al. 200r|). Dashed-dotted line in bottom panel corresponds to the [Mg/Fe]— a 
relation for Virgo cluster galaxies, by I Vazdekis et al.l d20oi) . Note that the slope of the [CN/Fe]— a relation increases with the density of 
the environment, but the [Mg/Fe]— er remains almost constant. 



density, the [Mg/Fe] ratio keeps constant and positive. We 
interprete the fact that [Mg/Fe] is positive and constant in 
terms of the formation timescales of the stars in ellipticals, 
as a function of the environment. 

If we consider the formation timescales for Mg and Fe 
stated above we can estimate that, whatever the environ- 
ment is, all stellar populations in ellipticals should be al- 
ready formed in less than ~ 1 Gyr, since this is the time 
needed by SNe la to pollute the ISM with iron-peak ele- 
ments. But, if we consider species with less separated for- 
mation timescales, such as Fe and CN, we do find differences 
in the ratio [CN/Fe] as a function of the environment, for 
a given a. That suggests a different formation timescale for 
the stars in ellipticals related to the environment proper- 
ties, in the sense that the more dense the environment is the 
shorter the formation timescales of the stars are (with an up- 
per limit of ~ 1 Gyr). This is in agreeement with a number 
of recent results (e.g. | Carretero et al.l |2004 iThomas et al.l 
l2005l ; lBernardi et al.ll2006h . 

Second, for our sample of galaxies in very rich clus- 
ters, we do find positive trends between [CN/H], [Mg/H], 
[CN/Fe], [Mg/Fe] and a. On the contrary, no correla- 
tions are found between [Fe/H] or Log(age) with the 
velocity dispersion. These relations are explained within 
a frame work where more massive galaxies produce more 
metals (jKodama fe Arimotol 1 19971 ; IChiosi fe Carrarol 120021 : 
Matteuccil 120031). O ur results are compatible with those of 
Trager et al,l (|2000r ) , who found a relation between the abun- 



dance of some a-elements, including Mg, and the velocity 
dispersion for a sam ple of old cluster ga l axies. Moreover, 
IVazdekis et al.l (|200ll ) and lYamada etafl (|2006f ) suggested 
that the CMR is driven essentially by metallicity and, in 
particular, by the abundance of the a-enhanced elements. 

The Log(age)— a relation presented in FigureUand the 
slopes of Figure [6jl show that the mean stellar age is nearly 
constant for all galaxies within rich clusters. This is not the 
case for the Virgo cluster galaxies, which present bigger ages 
for the most massive objects. This is due to the fact that 
smaller galaxies in low density envir onments present a larger 
scatt er in the age— a relation (see ISanchez-Blazquez et al.l 
2006). In fact, the slope in age— a for the Virgo cluster is 
also near zero if the sample is limited to galaxies with sigma 
larger than 200 km s _1 . 

Third, it is noteworthy that the slope of the relation 
[CN/Fe]— a increases with the cluster X-ray luminosity while 
the slope of the [Mg/Fe]— a is constant with the cluster 
mass (see Figure [6}. Because of that, we can only find 
environment-related [CN/Fe] differences in those galaxies 
with a < 300 km s _1 whereas, for more massive galax- 
ies, the [CN/Fe] values are almost equal independently of 
the cluster properties. This means that the ratio [CN/Fe] 
in most massive galaxies is less sensitive to the environ- 
ment than in intermediate- and low-mass galaxies, suggest- 
ing that star formation histories are less environment de- 
pendent as the galaxy mass increases. Therefore, [CN/Fe] 
appears to be an appropriate "chemical clock" for ellipticals 
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Figure 6. Slopes of the linear fits for the relations [CN/H]-ct (a), [Mg/H]-<r (b), [Fe/H]-tr (c), Log(age)-tr (d), [CN/Fe]-tr (e), 
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with the environment is found for the other elements. 



with a < 300 km s _1 but, when we move to very massive 
galaxies, the formation timescales differ so gently from one 
environment to another that this "clock" turns out to be 
insufficiently sensitive. In order to disentangle the different 
formation timescales between very massive galaxies in differ- 
ent environments, it is neccesary to find out a more accurate 
indicator. 



very short star formation histories, and which have been 
passively evolving since, at least, z ~ 2 . This is in agree- 
ment with a number recent results (e.g. lLabbe et alj|2005l ; 
iKriek et all l200fj| ; Ivan Dokkum fc van der Marell 120061 ). In 
comparison with less dense environments, the stars in ellip- 
tical galaxies of rich clusters must have formed at slightly 
earlier epochs and on a slightly shorter timescales. 



Our results suggest an overall picture in which massive 
elliptical galaxies in very rich clusters are old systems, with 
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5 CONCLUSIONS 

We have performed a detailed stellar population analysis of 
a sample of 27 massive elliptical galaxies (220 km s _1 < 
a < 450 km s _1 ) within 4 very rich Abell clusters: A115, 
A655, A963 and A2111. From the absorption indices H/3, 
CN2, Mg2 and Fe4383 we derived mean luminosity-weighted 
ages, relative abundances [CN/H], [Mg/H] and [Fe/H], and 
the abundance ratios [CN/Fe] and [Mg/Fe]. To obtain these 
results we used our new, high-resolution stellar population 
models (Vazdekis et al 2006, in preparation). We obtained 
the dependences of all parameters above on the galaxy ve- 
locity dispersion. 

Given the different formation timescales for each ele- 
ment, we have shown that the abundance ratios [CN/Fe] 
and [Mg/Fe] are the key "chemical clocks" for infering the 
star formation history timescales in ellipticals. In particu- 
lar, [Mg/Fe] provides an upper limit for those formation 
timescales, while [CN/Fe] apperars to be the most suitable 
"clock" to discriminate finely among them. This is true for 
those galaxies with a < 300 km s _1 as, for more massive 
galaxies, the [CN/Fe] ratio loses its "timing" sensitivity. 

Our [Mg/Fe] results are compatible with a scenario in 
which the stellar populations of elliptical galaxies, inde- 
pendently of their environment and mass, had formation 
timescales shorter than ~ 1 Gyr. This result implies that 
massive elliptical galaxies have evolved passively since, at 
least, 2 ~ 2. 

On the other hand, the information provided by the 
[CN/Fe] analysis suggests that, for a given galaxy mass (i.e. 
a given velocity dispersion), the star formation history must 
be shorter in those galaxies belonging to more dense environ- 
ments. Furthermore, we found that these timescales differ- 
ences become smaller for the most massive cluster galaxies 
(a > 300 km s -1 ), so they cannot be traced by the [CN/Fe] 
abundance ratio. 

Finally, we found that the mass-metallicity relation is 
steeper as the density of the environment increases. 

Our results do not allow us, at this point, to decide be- 
tween the two principal scenarios, hierarchical or monolithic, 
for galaxy formation, but do set certain strong constraints. 
Specifically, we have shown that galaxy assembly timescales 
for ellipticals are short and environment dependent, and this 
must be addressed by any valid theory of galaxy formation. 
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